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ABSTRACT 

Young Stellar Objects (YSOs) in the early evolutionary stages are very embedded, and thus they 
emit most of their energy at long wavelengths such as far- infrared (FIR) and submillimeter (Submm) . 
Therefore, the FIR observational data are very important to classify the accurate evolutionary stages 
of these embedded YSOs, and to better constrain their physical parameters in the dust continuum 
modeling. We selected 28 YSOs, which were detected in the AKARI Far-Infrared Surveyor (FIS), 
from the Spitzer c2d legacy YSO catalogs to test the effect of FIR fluxes on the classification of their 
evolutionary stages and on the constraining of envelope properties, internal luminosity, and UV strength 
of the Interstellar Radiation Field (ISRF). According to our test, one can mis-classify the evolutionary 
stages of YSOs, especially the very embedded ones if the FIR fluxes are not included. In addition, the 
total amount of heating of YSOs can be underestimated without the FIR observational data. 
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1. INTRODUCTION 

Stars are born in their parental molecular cloud 
cores through gravitational contraction. After the 
gravitational collapse begins, a protostar with its cir- 
cumstellar disk forms at the center of the collapsing 
core. Although the protostar/disk system is embedded 
in the dense infalling envelope at the early phase of 
evolution of star, the envelope gradually disappears as 
the protostar evolves to the main-sequence star. The 
Spectral Energy Distribution (SED) of a Young Stel- 
lar Object (YSO) is contributed by the protostar at 
Near-Infrared wavelengths whose temperature ranges 
from 2000 K to 4000 K, the disk at Mid-Infrared (MIR) 
wavelength which has a temperature gradient with ra- 
dius, and the envelope at Far-Infrared (FIR, 40 < A < 
350 /im) or Submillimeter (Submm, A > 350 /im) wave- 
lengths. Therefore, the analysis of the SED of a YSO 
can provide a way to identify the evolutionary stage 
of the YSO. In the SED analysis, however, we need to 
understand the effect of the interstellar radiation field 
(ISRF) since the envelope is heated by the ISRF as well 
as the internal source luminosity. 

There is a commonly known scheme classifying the 
evolutionary stages of YSOs based on the shape of the 
SED of a YSO (Lada 1987; Andre et al. 1993). In 
this proposed scheme, the Class object is in the earli- 
est stage after the collapse begins. The Class object 
already harbors a central heating source, but most of 
the mass is still in the dense envelope of gas and dust. 
Therefore, all emission from the central source is repro- 



cessed through the envelope, which is almost isothermal 
with T-30 K (Young et al. 2003). 

Class I object is still heavily embedded, but emission 
from the central protostar and disk escapes through the 
outflow cavities. As a result, the peak of the SED shows 
at MIR wavelengths. 

As the envelope material accretes to the central star 
through the disk and/or is cleared out by the outflows, 
the protostar/disk system becomes the dominant mass 
reservoir. Therefore, Class II objects are less embedded 
and have the emission excesses at near and mid-infrared 
wavelengths by dusty disk. The Class III objects emit 
energy mainly from central protostars and have only 
a small amount of gas in disks. The Class II and the 
Class III are often classified as classical T Tauri Stars 
and weak-line T Tauri stars, respectively (White et al. 
2007). 

As mentioned, in the early embedded phases, the ra- 
diation from central source undergoes the absorption/re- 
emission processes by the heavy envelope. As a re- 
sult, significant energy is emitted at FIR wavelengths. 
Andre et al. (2000) noted that deeply embedded Class 
0, whose energy peak is located at beyond 100 /^m, were 
often missed by IRAS observations unlike the Class I 
objects, which reveal the peak energy emission at the 
wavelength from 12 to 100 /xm. Therefore, the FIR data 
are crucial to study the evolutionary stages of embed- 
ded sources. 

The envelope is also heated by the ISRF so that the 
SED at FIR and Submm is affected by the energy in- 
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put from the ISRF as well as the central source. Since 
the heating from the ISRF also affects the total lumi- 
nosity of a YSO and the shape of the dust temperature 
profile, it is very important to constrain the strength 
of the ISRF, and thus the internal luminosity of the 
central heating soiircc. The external heating by the 
ISRF compared to the internal heating is not negligi- 
ble especially in the early embedded phase and in low 
luminosity sources (Dunham et al. 2008). 

The AKARI/Far-Infrared Surveyor (FIS) All-Sky 
Survey Bright Source Catalog combined with prece- 
dented catalogs such as those by the Spitzer obser- 
vations provides a useful tool to investigate the im- 
portance of the FIR fluxes in the calculations of the 
evolutionary indicators of YSOs, such as Tboi and 
Lboi/Lsmm- In this paper, AKARI/FIS All Sky Sur- 
vey Bright Source Catalog is used to investigate the 
importance of the FIR fluxes on the evolutionary in- 
dicators of YSOs, Tjjoi and Lboi/Lgmm- In addition, 
we model the SED of one YSO to study how the FIR 
fluxes can constrain the internal luminosity and the 
ISRF strength around the YSO. 

This paper is structured as follows. In section 2 
and 3, the catalogs used in this study and the crite- 
ria adopted to select our sample YSOs are introduced. 
The evolutionary indicators of YSOs are summarized 
in section 4. We calculate those evolutionary indica- 
tors of our samples and present the results in section 5. 
The dust continuum modeling of the SED of one sam- 
ple YSO is presented in section 6, and the summary 
and conclusions are followed in section 7. 

2. CATALOGS 

In this study, we use various infrared point source 
catalogs. All of them, except for one from the obser- 
vations with Spitzer Space Telescope, were obtained 
through the all sky survey. 

2.1 InfraRed Astronomy Satellite (IRAS) Cat- 
alog 

IRAS was launched on January 26, 1983 with a 
cooled 57 cm telescope for all-sky survey in 4 bands 
(12 /xm, 25 ^m, 60 ^m, and 100 ^m). The typical 
aperture size in each band is 3.5', 3.5', 3.6', and 4.7', 
respectively (Wheelock et al. 1994). Through IRAS 
all-sky survey, 72% of the sky was observed with three 
or more hours-confirming scans, and 95% was observed 
with two or more hours-confirming scans (Neugebauer 
et al. 1984). 

2.2 Two Micron Sky Survey (2MASS) Catalog 

Observations of the Two Micron Sky Survey began 
in 1997 and were accomplished in 2001 with the two 
telescopes located in the northern and southern hemi- 
spheres. The survey was performed for the entire sky 
using photometry system of 3 bands, i.e., J (1.235 /xm). 



H (1.662 /im), Ks (2.159 /im). The aperture size in all 
bands is 2.5"(Skrutskie et al. 2006). 

2.3 The c2d Spitzer legacy YSO Catalogs 

The Spitzer Space Telescope (SST) was launched in 
2003 which has three instruments, i.e., the Infrared Ar- 
ray Camera (IRAC), the Infrared Spectrograph (IRS), 
and the Multiband Imaging Photometer for Spitzer 
(MIPS). The IRAC has four bands (IRACl, IRAC2, 
IRAC3, and IRAC4), which have central wavelengths 
at 3.6 yum, 4.5 /xm, 5.8 /itm, and 8.0 /xm, respectively, 
and the beam FWHM is 1.66", 1.72", 1.88", and 1.98" 
in each band. Two bands of the MIPS, MIPSl and 
MIPS2, have effective wavelengths at 24 ^m and 70 
yum, respectively, and the beam FWHM in each band 
is 6"and 18". The YSO catalogs used for this study 
were provided by the Spitzer Legacy Program "From 
Molecular Cores to Planet Forming Disks" (c2d; Evans 
et al. 2003). The c2d YSO catalogs provide the IRAC 
and MIPS fluxes of identifled YSOs in five nearby, 
large molecular clouds (Serpens, Lupus, Chamaeleon 
II, Ophiucus, and Perseus) and 58 small cores. 

2.4 AKARI (ASTRO-F)/FIS Bright Source 
Catalog 

AKARI, launched on February 2004, is a cooled 67 
cm telescope with two focal plane instruments : the 
Far Infrared Surveyor (FIS) and the Infrared Cam- 
era (IRC). AKARI/FIS observed around 94% of the 
whole sky twice or more with a higher spatial resolu- 
tion and a wider wavelength coverage than those of the 
previous all-sky mission, IRAS (Kawada et al. 2007). 
Due to the higher resolution of AKARI/FIS observa- 
tions, AKARI/FIS could resolve the background and 
the source more effectively than IRAS (Jeong et al. 
2007). The FIS performed the all-sky survey with 
4 photometric bands (N60, WIDE-S, WIDE-L, and 
N160) which have an effective wavelength at 65 fim, 
90 ^m, 140 ^m, and 160 /zm, respectively. The beam 
FWHM of each band is 37", 39", 58", and 61", respec- 
tively (Kawada et al. 2007). 

3. DATA CRITERIA 

The sample YSOs selected for this study satisiy 
three criteria: 1. They are listed in the c2d catalogs, 2. 
Their MIPS2 70 /um fluxes are greater than 1 Jy, and 
3. They are detected by the AKARI but not by the 
IRAS. 

We use YSOs from the c2d YSO catalogs as our 
study targets. Considering the sensitivity limit of 
the AKARI/FIS, only those YSOs with 70 /zm fluxes 
greater than 1 Jy are used. The Key catalog for this 
study is the AKARI/FIS Bright Source Catalog since 
the goal of our study is to understand how FIR fluxes 
affect the properties of YSOs derived from observa- 
tions. The IRAS catalog could be used to study the 
effect of the FIR, but we adopt the AKARI/FIS catalog 
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Table 1. 

Properties of YSOs selected for this study 



Index 


c2d Region^ 


Spitzer Source Name*^ 


xt.A.^^zUUU.U ) 


JJeci.^^zuuu.Uj 


AivArii iName 


1 


Ophiuchus 


J163133.4- 242737 


247.891 


-24.4594 


1631337-242734 


2 


Ophiuchus 


J163135. 6- 240129 


247.899 


-24.0247 


1631357-240129 


3 


Lupus I 


J153848.4-344038 


234.703 


-34.6759 


1538486-344033 


4 


Perseus 


J032538.8+304406 


51.4102 


30.7361 


0325384+304410 


5 


Perseus 


J032637.5+301528 


51.6560 


30.2584 


0326374+301530 


6 


Perseus 


J032834.5+310051 


52.1426 


31.0154 


0328342+310055 


7 


Perseus 


J032837.1+311331 


52.1533 


31.2267 


0328368+311336 


8 


Perseus 


J032845. 3+310542 


52.1875 


31.0966 


0328450+310548 


9 


Perseus 


J032923. 5+313330 


52.3481 


31.5575 


0329235+313327 


10 


Perseus 


J033121.0+304530 


52.8372 


30.7585 


0331209+304531 


11 


Perseus 


J033229.2+310241 


53.1213 


31.0448 


0332291+310241 


12 


Perseus 


J033312.8+312124 


53.3031 


31.3563 


0333127+312123 


13 


Perseus 


J033317.9+310932 


53.3258 


31.1597 


0333182+310935 


14 


Perseus 


J033327.3+310710 


53.3642 


31.1198 


0333274+310711 


15 


Perseus 


J034255.9+315842 


55.7321 


31.9785 


0342557+315843 


16 


Serpens 


J182815.3+000-243 


277.063 


-0.0444 


1828152-000240 


17 


Serpens 


J182858.1+001724 


277.242 


0.2896 


1828580+001723 


18 


Serpens 


J182906.2+003043 


277.277 


0.5116 


1829064+003042 


19 


Serpens 


J182932.0+011843 


277.384 


1.3119 


1829322+011843 


20 


DC275.9+1.9 


J094645.9-510610 


146.691 


-51.1028 


0946457-510610 


21 


DC302.1+7.4 


J124540.0-552522 


191.417 


-55.4212 


1245402-552516 


22 


LlOO 


J171603.2-205657 


259.013 


-20.9489 


1716030-205656 


23 


L1251 


J223031.9+751409 


337.633 


75.2351 


2230318+751406 


24 


L723 


J191753.7+191220 


289.474 


19.2054 


1917539+191219 


25 


Lupus III 


J160918.1-390453 


242.327 


-39.0816 


1609184-390454 


26 


Lupus IV 


J160044.5-415531 


240.185 


-41.9242 


1600443-415527 


27 


Lupus IV 


J160115. 5-415235 


240.315 


-41.8752 


1601157-415231 


28 


Lupus IV 


J160221. 6-414054 


240.589 


-41.6811 


1602213-414052 



^Name of the region where each source is located. 

^AU source names are preceded by the prefix "SSTc2d." 



Table 2. 

2MASS and Spitzer flux densities of YSOs selected for this study 



Index 


J(1.25 («m) 
mjy 


H(1.65 /»m) 
mJy 


Srs(2.17 f»m) 
mJy 


IRAC1(3.6 /»m) 
mJy 


IRAC2(4.5 /»m) 
mJy 


IRAC3(5.8 /l»m) 
mJy 


IRAC4(8.0 /»m) 
mJy 


MIPS1(24 /»m) 
mJy 


MIPS2(70 /l»m) 
mJy 


1 


323±6.84 


515±27 


6n2±13.3 


62n±34,8 


54n±29,7 


489±23,7 


673±33.9 


1690±158 


256n±246 


2 


1,00±0.05 


19.0±0.5 


137±2.65 


589 = 32,8 


957 = 52,7 


1250±62,7 


1450±81.8 


3230±301 


10400±983 


3 


2.15±0.2S 


3.84±0.6 


6,60±0,54 


3.95±0.40 


2.76±0.52 


9.58±0.75 


42.7±3.11 


114±11.2 


2920±279 


4 








3.05±0.22 


11.5±1.18 


16. Oil. 18 


20.4±1.60 


1790±167 


31000±2900 


5 


O.IS(-) 


0.49(-) 


l-19±0-n8 


3.62±0.23 


9.27±a.46 


11. lie. 54 


12.0±0.56 


396±36.7 


4300±406 


6 


1.59±a.09 


5.43±0.24 


10.5±0.31 


54.7±2.66 


115±6.86 


221±10.9 


242±12.2 


1160±109 


2700±254 


7 


0.58(-) 


1.18±a.l4 


lO.l(-) 


30.1±1.75 


89.2±5.47 


267±13,3 


722±41.3 


6860±1150 


53800±5020 


8 


0.06(-) 


0.45(-) 


0.72±0.10 


1.22±0.06 


2,73±0,13 


2,67±0,13 


2,69±0.13 


213±19,9 


1520±150 


9 








0.76±0.07 


2.26±0.13 


2.47±0.20 


3,03±0.15 


191±17,7 


1150±110 


10 








0.12±0.01 


0.85±0.05 


0.83±0.05 


0.62±0.04 


14.7±1.37 


3880±366 


11 


0.71±0.06 


2.69±0.15 


5.81±0.21 


11.1±0.66 


17.2±0.91 


23.3±1.13 


32.8±1.58 


181±16.7 


1130±109 


12 


22.4±0.39 


88.8±1.31 


167±2.77 


301±16.3 


435±26.9 


642±31.8 


1120±57.5 


4040±382 


3280±321 


13 








0.17±0.01 


6.40±0.32 


36.4±1.79 


102±5.14 


670±63.0 


11700±1190 


14 








2.43±0.12 


6.68±0.32 


7.93±0.37 


12.4±0.58 


1570±147 


5410±510 


15 


95.7±2.12 


132±3.52 


164±2.87 


398±23.0 


346±21.3 


375±18.4 


540±26.6 


168±15.7 


2430±2390 


16 


37.1(-) 


36.4±1.21 


95. 2(-) 


174±9.05 


177±9.28 


180±11,2 


191±11.8 


270±73.6 


1670±169 


17 


82.5±1.98 


120±3.41 


111±2.26 


52.5±2.97 


36.7±2.47 


31.2±1,68 


28,5±1,68 


9,74±0,92 


1040±101 


18 








8.05±0.41 


45.0±2.76 


93.9±4,77 


129±7.62 


1320±139 


7240±713 


19 


35.5±0.78 


123±3.62 


246±4.75 


481±24.5 


786±42.6 


1140±57,5 


1830±138 


4370±407 


5190±501 


20 








0.22±0.02 


0.94±0.06 


0.69±0.06 


0.59±0.04 


76.4±7.05 


6590±615 


21 








0.13±0.01 


1.10±0.05 


2.04±0.12 


2.39±0.11 


87.7±8.09 


1750±164 


22 


0.69±0.07 


3.88±0.30 


10.5±0.40 


86.7±4.66 


146±7.64 


203±10.7 


306±16.7 


1140±106 


2540±235 


23 








0.13±0.01 


0.42d=0.04 


0.31±0.04 


0.16±0.03 


4.97±0.48 


1400±131 


24 








0.37±0.02 


4,58±0,24 


11.5±0.55 


16.3±0.79 


295±27,3 


9040±838 


25 








0.25±0.03 


l,00i0,06 


0,98±0.07 


0,55±0,05 


32,4±3,04 


2610±252 


26 


263±6.29 


344±9.82 


305±6.19 


177±9.40 


141±7.04 


140±6.81 


213±11.2 


590±55,2 


1050±103 


27 


0.42±0.05 


1.67±0.09 


4.83±0.19 


8.36±0.44 


9.92±0.49 


8.98±0.43 


7.70±0.36 


75.9±7.03 


1220±123 


28 


1.63±0.10 


3.36±0.18 


4.66±0.21 


5.10±0.30 


4.87±0.25 


12.6±0.62 


42.1±2.03 


144±13.4 


1030±107 
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Table 3. 

AKARI FIS flux densities of YSOs selected for this study 













T 

Index 


nijy 


mJy 


mJy 


mJy 


1 


O O ^ < 1 /"I /~1 

3874±209 


3638±258 


(a) 




2 


12210±1140 


12420±1360 


14430±4470 


21650±5660 


3 


3261±505 


46G3±104 


5360±1720 




4 


80d80±10600 


133200±21000 


124d00±21200 


217400±28800 


5 


4453±593 


6092±420 


12880±2490 


10740±1110 


6 


3872±300 


3757±320 


6325±2010 


10490±947 


7 


101400±1490 


64350±11200 


58090±10400 


84570±885 


8 




2339±43 


6789±1820 




9 




1758±162 


7982±2490 




10 


3121±498 


9190±259 


20510±3760 


25870±3390 


11 




2431±275 




10070±182 


12 


3993±333 


4712±182 






13 


15260±1370 


25020±1530 


46290±4150 


57090±12500 


14 


4140±252 


A A A C\ 1 1 > 1 i'l 

4442 ±133 






15 


6654±754 


12640±1180 


34980±7450 


21750±3520 


16 


2298±158 


3590±138 


5456±1050 




17 




2574±99 


9914±1890 


7043±1550 


18 


38380±670 


41090±1970 


62110±12400 


109600±7900 


19 






23600±6020 


16030±4390 


20 


5544±652 


11470±452 


20810±2580 


19110±2430 


21 




2777±130 


9229±2700 




22 




4822±632 


14440±4450 


8270±1480 


23 




2164±101 


7670±530 




24 


8573±1530 


15800±582 


31160±2020 


34310±8050 


25 




3807±505 


7356±938 




26 




1470±236 






27 




2146±218 


6882±981 




28 




1542±66 







* Fluxes of the flux density flag smaUer than 2 are not considered. 



since the AKARI resolved the point sources effectively 
due to its higher resolution (Jeong et al. 2007) and 
absolute calibration errors arc expected to be less than 
20% in all bands (Yamamura et al. 2010). In addition, 
most of 100 yum fluxes for YSOs matched with IRAS 
do not have a good flux quality. Therefore, we have 
used only AKARI/FIS sources with a good flux qual- 
ity for our study. Besides, previous studies of those 
YSOs detected in AKARI but not detected in IRAS 
have been conducted without the FIR data. There- 
fore, it is meaningful to compare previous studies with 
our study, which includes FIR fluxes. 

Wc carry out cross correlation with 6" among cat- 
alogs based on the coordinates provided in the c2d 
YSOs catalogs considering positional uncertainties (Ya- 
mamura et al. 2010). There is little flux dependency of 
the position error on the current flux range of the sam- 
ple. The AKARI/FIS catalog provides the flux density 
quality flag, which indicates whether observed flux den- 
sity is reliable or not. For our study, we use only the 
fluxes flagged as "reliable" from "conflrmed" sources 



(Yamamtira et al. 2010). Based on these criteria, 28 
YSOs have been selected as our study sample, as listed 
in Table 1. The flux densities and their errors of these 
28 YSOs are presented in Table 2 and 3. 

4. EVOLUTIONARY INDICATORS 

In order to classify the evolutionary stages of YSOs, 
a few evolutionary indicators have been used. We sum- 
marize them below. 

4.1 Spectral Index 

Since Lada & Wilking (1984) recognized that the 
SEDs of YSOs were falling into natural groups, Lada 
(1987) first codified the tripartite class system using 
the spectral index (a) between 2 and 20 /xm: 



a 



rflog(A) 



(1) 
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Table 4. 

Physical properties of the selected 28 YSOs 



Index 


Spectral Index^a) 


T ( T ^ 1 

L'int'y'-'Q)) 




{R} 


Lboi{LQ) 


Without 
c IH riuxes 


with 
b IK nuxes 


without 
FIR fluxes^ 


with 
FIR fluxes^ 


1 


-0.610 


3.371 


1202 


1143 


0.763 


0.804 


2 


0.140 


12.590 


452 


380 


0.919 


1.116 


3 


0.580 


3.814 


180 


119 


0.110 


0.199 


4 


2.160 


35.148 


76 


53 


3.094 


12.593 


5 


1.090 


5.488 


105 


75 


0.490 


0.894 


6 


0.830 


3.544 


307 


234 


0.846 


1.167 


7 


2.350 


59.015 


119 


85 


6.997 


13.863 


8 


1.110 


2.065 


114 


84 


0.196 


0.282 


9 


1.510 


1.588 


104 


67 


0.158 


0.360 


10 


0.980 


4.983 


56 


41 


0.316 


0.856 


11 


0.400 


1.562 


283 


145 


0.191 


0.435 


12 


0.410 


4.255 


496 


458 


2.822 


3.078 


13 


3.330 


14.064 


92 


61 


1.211 


2.973 


14 


1.930 


6.811 


114 


107 


0.962 


1.066 


1 l^ 
10 


-u.you 




±u±o 


o^y 


1.300 


2.489 


16 


-0.850 


2.256 


847 


641 


0.730 


0.981 


17 


-1.920 


1.445 


1559 


1033 


0.466 


0.711 


18 


1.700 


8.957 


138 


64 


1.199 


6.191 


19 


0.260 


6.550 


546 


448 


4.168 


5.168 


20 


1.660 


8.199 


59 


47 


0.798 


1.640 


21 


1.910 


2.357 


76 


50 


0.246 


0.502 


22 


0.670 


3.346 


334 


240 


0.710 


1.037 


23 


0.700 


1.911 


56 


36 


0.164 


0.311 


24 


2.150 


11.036 


72 


53 


1.200 


2.552 


25 


1.100 


3.432 


60 


47 


0.141 


0.227 


26 


-0.620 


1.458 


1492 


1430 


0.491 


0.513 


27 


-0.190 


1.679 


232 


136 


0.051 


0.098 


28 


0.920 


1.432 


276 


323 


0.059 


0.048 



^ The iiitcnal luminosity calculated from the flux at 70 /nm (Dunham at al. 2008). 
^ Tbol a^nd -^'boi calculated only with the c2d data. 

^ Tiai and L^oi calculated with the c2d data combined with the AKARI/FIS data. 



where A is the wavelength and S\ is the flux density 
at that wavelength. The boundaries of class system by 
Lada (1987) are as follows: 

I 0<a^3, 

II -2 < a ^ 0, and 

III -3 < a ^ -2. 

After the tripartite class system by Lada (1987) was 
introduced, the class system has been improved by 
Wilking et al. (1989), Andre et al. (1993, 2000) and 
Greene et al. (1994) who formalized the 4-class system, 
as follows: 

I 0.3^ a. 

Flat -0.3 ^ a < 0.3, 

II -1.6 ^ a < -0.3, and 

III a < -1.6. 



Although this class system was based on the slope 
between 2.2 and 10 /xm, they asserted that these classes 
showed no systematic deviations from classes based on 
the slopes between 2.2 and 20 /im. The c2d catalogs 
provide a least-squares fit to all photometry between 
2 and 24 /im to determine the spectral index {a) and 
classify objects into the four classes defined by Greene 
et al. (1994) (Evans et al. 2009). We use a provided 
by the c2d catalogs as an evolutionary indicator of our 
selected YSOs. 

4.2 Bolometric Temperature 

Myers & Ladd (1993) suggests T^oi as a measure of 
the circumstellar obscuration and evolutionary devel- 
opment of a YSO. We calculate the bolometric temper- 
ature {Tbol), which is the temperature of a blackbody 
having the same mean frequency P as the observed con- 
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Fig. 1. — The plot of a vs. Ttoi for a subset of the sources 
from Table 4. We present sources only with Tboi < 200 K 
in this figure. The numbers in this figure indicate the in- 
dex numbers of YSOs designated in Table 1, 2, and 3. The 
horizontal dashed line indicates the boundary between the 
embedded (class O/I) and less embedded (class II) sources, 
and the vertical dashed line indicates the boundary be- 
tween class {Tboi < 70 K) and class I (70 K < Thoi < 
650 K). The left panel shows the values of Tboi calculated 
without the AKARI FIR data, and the right panel shows 
those calculated with the AKARI FIR data as well. The 
numbers of red color denote the YSOs whose evolutionary 
stage changed from class I to class in this comparison. 



tinuum spectrum, using the prescription given in Myers 
& Ladd (1993), 

Tboi = [C(4)/4C(5)]/ii?/fc = 1.25 x 10-^^i?KHz-\ (2) 

where C (n) is the Riemann zeta function of argu- 
ment n, h is Planck's constant, and k is Boltzmann's 
constant. The mean frequency 9 is the ratio of the first 
and zeroth frequency moments of the spectrum. 



(3) 



Chen et al. (1995) showed that evolutionary classes, 
including class 0, ranged depending on T^oi as follows: 

Tboi < 70, 

1 70 ^ Tboi ^ 650, and 
II 650 < Tboi <, 2800. 

4.3 The ratio between Lboi and Lgmm 

The ratio of bolometric to submillimeter luminosity 
{Lboi/ Lsmm) can be an evolutionary indicator of YSOs, 
especially in the early embedded phase (Andre et al. 
1993). Lboi is calculated by integrating over the full 
observed SED, 



Lboi = / 47rL) Sudv, 



(4) 



while the submillimeter luminosity is calculated by in- 
tegrating over the observed SED for A > 350 ^m, 



iy=c/350/^m 



AirD'^S.Aiy. 



(5) 



where Si, is the flux density and D is the distance. 
Young & Evans (2005) suggested these boundaries 
for classifying protostars; Lboi/ Lgmm — 35 for Pre- 
Protostellar Core/Class 0, Lboi/ Lsmm = 175 for Class 
O/I, and Lboi/Lsmm ^ 2000 for the Class I/II transi- 
tion. 

The spectral index a is calculated without the FIR 
data so that it cannot be used to subdivide the em- 
bedded stage into Class and Class I. However, Tboi 
and Lboi / Lsmm have ability to distinguish Class from 
Class I. A drawback to using Lboi/ Lsmm is in the lack 
of the complete data set at A > 350 fim in many YSOs. 

5. RESULTS 

We calculated Tboi and Lboi of 28 YSOs selected by 
the criteria presented in section 3 in both cases of (1) 
only with NIR and MIR and (2) including the AKARI- 
FIR data as well. Table 4 lists the spectral index pro- 
vided from the c2d catalogs and Tboi and Lboi calcu- 
lated by integrating the observed flux density within 
the wavelengths covered by observations. 

We calculated the internal luminosity (Lint) listed 
in Table 4 using the equation presented in Dunham et 
al. (2008): 



L™t =3.3 X lO^^^^V^L©, 



(6) 



where is a flux density at 70 /im in cgs units (ergs 
cm~^ s~^) and normalized to 140 pc. This equation is 
derived using the results of the linear least-square fit to 
the modeled Lint versus f 70 in log-log space. Since F70 
is an observable quantity, the correlation they found 
between F70 and Lint can be used to estimate Lint for 
sources which either lack sufficient data to constrain 
radiative transfer models or have such data but have 
not yet been modeled. 

We compare the values of Tboi of the 28 YSOs calcu- 
lated with and without the AKARI FIR fluxes in Fig. 1. 
According to a, 20 out of 28 YSOs are classified as 
Class O/I. As mentioned, in the previous section, how- 
ever, a cannot subdivide the embedded phase. There- 
fore, for the detail study of the early embedded YSOs, 
we need to use Tboi and Lboi/ Lgmm- As presented in 
Fig. 1, class and class I objects can be divided based 
on Tboi] Tboi < 70 K for class and Tboi > 70 K for 
class I. When the FIR fiuxes are not included in the 
calculation of Tboi (the left panel of Fig. 1, see Table 4 
for the whole samples), 5, 19, and 4 YSOs are classified 
to class II, class I, and class 0, respectively. However, 
when we include the FIR fiuxes in the calculation of 
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Fig. 2. — The plot of Lboi vs. Tboi for the sources from 
Table 4. The numbers denote the index of YSOs as listed 
in Table 1. The categories of two panels are the same as in 
Fig. 1, that is, without (left panel) and with (right panel) 
the FIR fluxes in the calculation of Lboi- The dashed line 
indicates the boundary between class and class I in Ttoi- 
The numbers of red color denote the YSOs whose Ltoi in- 
creases more than twice when the FIR fluxes are included 
in the calculation of Ltoi ■ 



Ttoi (the right panel of Fig. 1), 3, 15, and 10 YSOs are 
classified to class II, I, and 0, respectively. 

That is, 6 YSOs are reclassified to class from class 
I, and 2 YSOs are reclassified to class I from class II. 
As a result, YSOs detected in the FIR wavelengths 
are mostly (25 out of 28) embedded. This compari- 
son shows that the FIR fluxes are very important to 
know the evolutionary stages of YSOs, especially for 
the embedded ones. 

Fig. 2 and 3 present the same comparison but for 
Lboi ■ When the FIR fluxes are included in the calcula- 
tion of Lboi, the values of Ltoi of 10 YSOs increase by 
more than a factor of 2. Two of them (sources 4 and 
18) even increase by more than a factor of 4 in their 
values of Lboi- This change in Lboi possibly increases 
the ratio, Lboi/Lgmm, resulting in different evolution- 
ary stages of the YSOs. According to our comparisons 
in the calculations of Tboi and Lboi with and without 
the FIR fluxes, we can see how important knowing the 
FIR fluxes of YSOs, especially the embedded ones, is. 

6. ONE-DIMENSIONAL DUST RADIATIVE 
TRANSFER MODELING 

The original radiation from the central source and 
circumstellar disk is scattered, absorbed and reemitted 
by the dusty material in the envelope. The embed- 
ded sources have thick envelopes, where most of mass 
resides, and which can be well described by a sphere. 
Therefore, in order to get the detailed information on 
YSOs, the emerged SED must be modeled. In this 
study, we modeled the SED of a YSO using a 1-D ra- 
diative transfer code to test how the FIR fluxes can 
constrain the physical parameters of the YSO. 



Fig. 3. — The same plot as Fig. 2 but for YSOs with 
higher Lboi- The values of Lboi of sources 4 and 18 increase 
more than 4 times when the FIR fluxes are included in the 
calculation of Ltoi- 



The selected YSO for this test is the source 13 in 
Table 1. Its Spitzer source name is J033317.9-I-310932, 
and it is located in a large molecular cloud, Perseus. As 
listed in Table 4, Lboi is '^1.2 and ~3.0 Lq when cal- 
culated without and with the FIR fluxes, respectively. 
In the early embedded phase, the heating by the ISRF 
cannot be ignored compared to the heating by the cen- 
tral source (Young et al. 2003). Therefore, in order to 
know the accurate internal luminosity, which is a cru- 
cial quantity to understand the formation mechanism 
of the protostar, separately from the ISRF luminosity, 
we must model the observed SED. 

J033317.9-h310932 was not detected by the 2MASS 
observations. This source is listed as Bolo 80 in the 
Bolocam survey of Perseus for 1.1 mm dust continuum 
emission (Enoch et al. 2006). The observed fluxes and 
calculated Tboi and Lboi of this YSO are listed in Ta- 
ble 2, 3, and 4. We use the one-dimensional radiative 
transfer code, DUSTY (Ivezic et al. 1999) to model the 
SED of J033317.9+310932. The program, ObsSphere 
(Shirley et al. 2002) is then used to match the results 
from DUSTY with the actual observed fluxes at wave- 
lengths longward of 100 /im. ObsSphere calculates the 
fluxes emerging only from the envelope, which is the 
dominant flux source at long wavelengths, based on 
the physical conditions derived from DUSTY. Refer to 
Dunham et al. (2006) for the detailed description of 
the model input parameters. 

We modeled the SED with and without the FIR 
fluxes and compared the derived properties of the cen- 
tral heating source (star/disk) and the envelope. The 
envelope is assumed to have a power-law radial den- 
sity profile [n{r)oor^P]- We have tested the envelope 
model by changing the p parameter from 1.5 to 2.0. 
According to the inside-out collapse model (Shu et al. 
1977), the density profile has the power of 1.5 and 2 at 
radii smaller and larger than the infall radius, respec- 
tively. The heating sources for our model are the in- 
ternal star/disk system and the ISRF. In general, since 
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Fig. 4. — The best-fit model of tfie SED of 
J033317. 9+310932 when only the Spitzer data are used. 
The green solid line indicates the best-fit model SED ob- 
tained by DUSTY, and red open diamonds represent the 
Spitzer photometry data. The blue dashed line shows the 
input SED for the model. The purple crosses show the 
emerging SED only from the envelope which is calculated 
with ObsSphere based on the model density profile and the 
temperature profile obtained from DUSTY. The star sym- 
bols denote the AKARI fluxes, which are not included for 
this SED model. The best-fit model only with the Spitzer 
fluxes does not match well the AKARI fluxes. 



cores with YSOs are embedded in bigger clouds, the 
ISRF reaching the surface of a core is attenuated by 
the surrounding material. We tested various attenua- 
tor of the ISRF, which is parameterized by the visual 
extinction (Ay) since the UV field strength of the ISRF 
relative to the local ISRF, Go is given by e"^ ®^'' (Lee 
et al. 2004). 

For the dust opacity, as many studies did (Shirley 
et al. 2002; Young et al. 2003; Evans et al. 2005), we 
adopted the 0II5 dust model, which is for the thin ice 
mantles coagulated for 10^ years at a gas density of 10^ 
cm~^ (Ossenkopf & Henning 1994). 

To find the best-fit envelope model we calculated the 
reduced (Dunham et al. 2006), 

a,(A.)2 ' 

where S°^^{\i) is the observed flux density at Ai, 
S^°'^{Xi) is the simulated flux density from this model 
at same wavelength, and a^{Xi) is the uncertainty in 
the observed flux density. And n is the number of the 
data points, so the model with the FIR fiuxes has four 
more data points than the other model without the FIR 
fluxes for the case of J033317.9-h310932. 

Fig. 4 shows the best-fit model without the FIR 
fluxes. The model parameters of this model are listed 
in Table 5. The best-fit model does not match very 



Fig. 5. — The best-fit model SED when the AKARI FIR 
fluxes as well as the Spitzer data are included. The purple 
triangle and open blue diamond indicate the observed 1.1 
mm Bolocam flux and the flux simulated from the model 
with the aperture used for the 1.1 mm flux measurement, 
respectively. Other colors and symbols are the same as in 
Fig. 4. 



well the observed 8 and 24 /im fluxes. This is possi- 
bly because the disk emission is not well modeled with 
our 1-D calculation. As seen in Fig. 4, this bcst-flt 
SED cannot match the AKARI fiuxes. According to 
this best-fit model, the ISRF is completely obscured 
by surrounding medium, and the internal luminosity is 
about 3.8 Lq. 

Fig. 5 presents the best-fit model when the FIR 
fiuxes are included. The parameters of this model are 
also listed in Table 5. The internal luminosity increases 
~25% compared to the model without the FIR fiuxes, 
and a higher UV field strength of the ISRF is neces- 
sary to fit the AKARI fluxes. The model SED shows 
that the peak of the SED is located at the wavelength 
of ^160 /im, which is the longest wavelength of the 
AKARI FIS. In order to constrain the SED better, 
however, observations at longer wavelengths than the 
AKARI bands are necessary. 

The density (upper panel) and temperature (lower 
panel) distributions of the envelope with radius for the 
best-fit models with (blue solid line) and without (red 
dashed line) the FIR fiuxes are presented in Fig. 6. 
Since AKARI FIR fluxes listed in the catalogue are 
not color corrected ones, we applied color correction to 
AKARI FIR fiuxes by using our modeled SEDs. The 
variation of color correction factor in our models does 
not exceed 10%. The envelope density profile of the 
best-fit model with the FIR fiuxes has the power of 
1.9, which is similar to the density profile of the static 
envelope in the Shu's inside-out collapse model. This 
result might indicate that J033317. 9+310932 is in a 
very early stage of collapse, so the infall radius is still 
small. 

The physical parameters for the best-fit model with 
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Table 5. 

The best-fit models for the SED of J033317.9+310932 



model 




R 2 

dinner 


R 3 


M ■* 






r 7 


rp 8 

star 






AU 


AU 


Mq 






Lq 


K 


without FIR fluxes 


1.6 


650 


20000 


5.9 


1.5x10*^ 





3.8 


2800 


with FIR fluxes 


1.9 


880 


29000 


7.6 


2.4x10'' 


0.4 


4.7 


3000 



^The power of the power law envelope density profile. 
^The inner radius of the envelope. 
^The outer radius of the envelope. 
*The mass of the envelope. 

5 The fiducial density at 1000 AU. 

6 The UV field strength of the ISRF relative to the local ISRF. 

^The internal luminosity caused by the accretion to the star/disk system. 
*The temperature of the central protostar. 
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Fig. 6. — The envelope physical structures obtained from 
the dust continuum modeling as a function of radius. The 
upper and lower panels present the density and temperature 
structures, respectively. The red dashed lines are for the 
best-flt model without the FIR fluxes while the blue solid 
lines are for the best-flt model with the FIR fluxes. 



the FIR fluxes show differences from those for the best- 
fit model without the FIR fluxes. As mentioned above, 



the power of the density profile is different in the two 
models. In addition, the envelope mass of the best- 
fit model with the FIR fluxes is about two solar mass 
greater than the best-fit model without the FIR fluxes. 
Other important parameters are the internal luminos- 
ity (Lint) and the UV strength of the ISRF (Gq), which 
are two heating sources of the envelope. According to 
our test. Lint and Go, constrained by the SED model- 
ing, are probably close to L^t = 4.7 Lq and Gq = 0.4 
rather than Lint = 3.8 Lq and Go = 0. Here, Go = 
0.4 is equivalent with the UV field attenuated by Ay 
= 0.5 mag. 

7. SUMMARY AND CONCLUSIONS 

We tested the effect of FIR fluxes in classiflcation 
of the evolutionary stages of YSOs and in calculation 
of the internal luminosity as well as the UV strength 
of the external heating source, the ISRF. Twenty eight 
sample YSOs, which have their counterparts only in the 
AKARI/FIS Bright Source Catalog, were selected from 
the c2d YSO catalogs. First, we calculated an evolu- 
tionary indicator, T},oi, and found that 6 objects, which 
are classifled into class I when only Spitzer fluxes were 
used, were reclassified into class when the AKARI 
FIR fluxes were included. Two objects, which were 
classified into class II when only the Spitzer data were 
included, were also reclassified into class I. This result 
shows that the FIR observations are very important in 
the classification of the evolutionary stages of the early 
embedded YSOs since they emit most of their energy 
at the long wavelength regime. 

In order to test how the FIR fluxes constrain the 
physical structure of the envelope, the internal lumi- 
nosity, and the UV strength of the ISRF, we modeled 
the SED of one of our sources. In the early evolution- 
ary phase, the external heating by the ISRF can not 
be ignored compared to the internal heating. 
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According to our dust continuum modeling, when 
the FIR fluxes are included, the density profile of the 
envelope is steeper (p=1.9 compared to p=1.6) indicat- 
ing most of envelope material is possibly static, that 
is, the infall radius is very small in the sense of Shu's 
insidc-out collapse model. The estimated internal lu- 
minosity is constrained to 4.7 Lq when the FIR fluxes 
are included (compared to 3.8 Lq obtained when the 
FIR fluxes are not included). The UV strength of the 
ISRF is Go=0.4 (attenuated by Ay = 0.5 mag) in the 
model with the FIR fluxes, but in the model without 
the FIR fluxes, the ISRF needs to be completely atten- 
uated (Go=0). Our modeling result also indicates that 
fluxes at FIR and even longer wavelengths are crucial 
to better constrain the physical parameters of embed- 
ded YSOs. 
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